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Abstract 

In the field theoretical Glauber-Gribov framework in the eikonal approxi- 
mation, we show that, in hadron-nucleus collisions in the high-energy limit and 
under specific conditions, the rescattering contributions to the J/ip produc- 
tion cross-section vanish, leading to an A 1 dependence on the atomic number 
A - except for nuclear effects in the nucleus structure function, the so-called 
shadowing. We show that the RHIC data on J ftp production in dAu collisions 
can be described with nuclear shadowing alone (cr a bs = 0). The corresponding 
results for Au Au collisions are also presented and compared with available 
data. 
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1 Introduction 



A surprising result in J ftp production in hadron-nucleus (or deuteron-nucleus) col- 
lisions is the substantial decrease of the so-called nuclear absorptive cross-section, 
Cabs , between CERN and RHIC energies [1, 2]. In this paper we propose a scenario 
which explains the observed decrease of a a b s with increasing s. 

We are going to proceed as follows: In Section 2 we develop our formalism for 
the absorptive cross-section a a b s in the framework of Glauber-Gribov field theory. 
In Section 3 we introduce a dynamical, non linear nuclear shadowing. It is a con- 
sequence of multiple scattering, and it is determined in terms of diffractive cross 
sections. In Section 4 we present our results for dAu and AuAu collisions at RHIC 
energies. We will finish by our conclusions in Section 5. 

2 The absorptive cross-section 

We proceed in the framework of Glauber-Gribov [3] field theory in the eikonal ap- 
proximation. We shall consider rescattering of both the light partons of the projec- 
tile, with a typical hadronic cross-section a, and the heavy cc system, with cross- 
section a. We shall assume that the dominant contribution for the J/ip production 
is the conventional mechanism such as gluon-gluon fusion. In the field theoretical 
language, we assume that the J/vfj is produced in a cut interaction of the projectile 
with a target nucleon. Production in the blob of a diagram is assumed to be small. 

For simplicity we shall assume throughout this paper that the amplitudes of both 
the projectile and heavy system scattering are purely imaginary. Let us denote by 
<tX(c) the pA inelastic cross-section, at fixed impact parameter, calculated in the 
Glauber model with a pN = a, i.e. 1 

= 1 - [1 - ° T A {b)] A = £ (fj^nib) [1 - vT A (b)] A - k . (1) 

Here T^(b) is the nuclear profile function normalized to unity, and k is the number 

of cut interactions. Since we allow interactions of both the projectile, with cross- 

1 For simplicity we have used the notation cross-section for this quantity. In fact it corresponds 
to a dimensionless quantity, the differentiated cross-section da/d 2 b. This is also valid for the later 
equations, where we have used this kind of notation too. 
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section a, and the heavy system, with cross-section a, to take place, the cross- 
section of the J/ip production is proportional to o l A {o + a) [4]. However, since we 
require that at least one interaction is cut, only its inelastic part o"X(°" + ^0 does 
contribute. Moreover, in our mechanism the J / ip is produced in a cut interaction of 
the projectile. Thus, we have to substract from a™ (a + a) the contribution with no 
cut interaction of the projectile, i.e. 

E [fj [(«■ + °) k - **] T » (1 - (o- + Z)T A {b)) A - k = . (2) 

Without the second term in the squared bracket (proportional to a k ) we have just 
the expression of a™ (a + a). This second term has been substracted since at least 
one cut interaction of the projectile (in which the heavy system is produced) has to 
be present. We see from (2) that this difference is equal to cr™(b). This means that 
all eikonal type rescatterings of the heavy system have cancelled out. 

The origin of this result is a cancellation between different discontinuities (cut- 
tings) of the diagrams in the Glauber-Gribov theory which is contained in the eikonal 
model - but has a more general validity. 

Let us now decompose a into two terms, a small one corresponding to the pro- 
duction of J/ip and the rest in which the J/ip is not produced : a = a^ N + '. 
From eq. (2), it is clear that the contribution to the J ftp cross-section can be written 
as 

= E [fj [° k ~ KW] [1 - T A (b)a)] A - k = a%(at N ) - A^ N T A (b) . 

(3) 

We see that the same cancellation is at work. We see from eq. (3) that only rescat- 
terings of eikonal type involving the small piece, a^ N , of a are left (self absorption) 
[5]. The conventional shadowing, related to diffractive cross-section <td, is due to 
diagrams of a non-eikonal type involving the triple Pomeron coupling and it is not 
included in the above formulae. It can be easily incorporated, as we will do below. 

Apart from these conventional shadowing corrections, we see from eq. (3) that 
J/ijj production is proportional to A 1 . We shall show below that this result is valid 
only at high energy. It is in sharp contrast with the situation at low energies where a 
rescattering of the pre-resonant cc system with the target nucleons takes place, with 
cross section a abs = (1 —e)a. a abs is the cross-section of the interaction of the cc pair 



with a target nucleoli in which the cc is transformed in such a way that it has no 
projection into J/ip (open charm is produced instead). We can write a a b s = (1 — e)a 
where a is the total cc — N cross-section and ea (e < 1) is the contribution to a of all 
intermediate states having a non-zero projection into J/ip. Since this contribution 
is expected to be quite small we should have e <C 1 for J/tp and, on the contrary, 
e ~ 1 for open charm production. Thus, for open charm a abs = (1 — e)a ~ while 
for J/ip production <r a 5 S ~ a. 

The origin of the change in the physical picture between low and high energies 
is well known [6] [7] [4]. The mass M c ^ of the heavy system introduces a new scale, 
the coherence length l c 

1 M- 

A = - = m p ^ . (4) 

l c S Xi 

At high energies the coherence length is large and the projectile interacts with 
the nucleus as a whole. At low energies the coherence length is small and the 
probabilistic picture with a longitudinal ordering in z of the various interactions is 
valid. In the first interaction at z the heavy system is produced and in successive 
interactions, at larger values of z, it rescatters with nucleons along its path. The 
change with energy of the rescattering mechanism has a clear physical interpretation. 
Due to the presence of the heavy system some contributions to the production cross- 
section (corresponding to some particular discontinuity of a given diagram) contain 
a non-vanishing minimal momentum transfer (t min ^ 0) and are suppressed by the 
nuclear form factor. 

Following ref. [4] the effects of the coherence length can be taken into account 
by writing the total pA cross-section as 

A n 



tot 



(i) = EL ( 5 ) 



n=l \ U / j=l 



where n is the number of interactions (both cut and uncut) of both the projectile and 
the heavy system. Here we have to consider a particular ordering of the longitudinal 
coordinates Zi of the n interacting nucleons : z± < z 2 < • • • < z n , and the n-th power 
of the nuclear profile function T^(b) is replaced by 

r+oo f+oo r+oo n 

T^(b) =n\ dzj dz 2 --- dz n cos(A(^ - z 3 )) J] (6) 

J -co Jzi Jz n -i i=l 
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where A is the inverse of the coherence length, eq. (4), and pa is the nuclear density. 
Note that for A = 0, corresponding to asymptotic energies, all integrals T^\b) are 
equal to T^{b). For A ^ 0, the value of j in eq. (5) depends on the discontinuity of 
the diagrams we are considering. Namely, interactions with nucleons 1 to j — 1 are 
uncut and located to the left of the cutting line. Interaction j can be either cut or 
located to the right of the cutting line. All other interactions from j + 1 to n can be 
cut or located to the right or to the left of the cutting line - all possibilities have to 
be added. To this contribution we have to add its complex conjugate, which gives 
just a factor of 2 for purely imaginary amplitudes. 

It is now straightforward to obtain a^p in eq. (5). We have to remember that 
when an interaction is uncut (located to either the right or the left of the cutting 
line), its amplitude (ia = —a/2, for purely imaginary amplitudes) is unchanged. 
When the cutting takes place through the interaction, its amplitude is replaced by 
ea in the case of an interaction of the heavy system. Indeed, only the fraction ea 
of the cc — N cross-section, corresponding to the intermediate states with non-zero 
projection into J/ip, does contribute 2 . In the case of an interaction of the projectile 
its amplitude is replaced by a. However, in order to keep track of the cut interactions 
of the projectile, we shall replace it by 5a with 8—1. 

Proceeding in this way it is easy to obtain the expressions for a^ in eq. (6) [4]. 
We get for j > 1 

= 2 (-| - \ ) [| + \ + v{6 - 1) + He - 1)] (a(8 - 1) + a(e - l))""' 

(7) 

and for j = 1 

aV = (a6 + ve)(v(6-l) + <T(e-l)r 1 . (8) 

The first factor of eq. (7) corresponds to the j — 1 uncut interactions. The second 
one corresponds to the jth interaction which is either cut (which gives 8a for the 
interaction of the projectile and ea for that of the heavy system) or located to the 
right of the cutting line (which gives — (a + a)/2). The last factor corresponds to 
the interactions j + 1 to n. For j = 1 the formula is different. Since a diagram can 



2 The introduction of e is most important at low energies where the heavy system undergoes 
successive interactions with nucleons in its path and has to survive all of them. It does not appear 
in the formula for s — > oo. 



5 



only be cut either through or in between interactions, the first interaction has to be 
cut and there is no complex conjugate term - because all other interactions have to 
be located to the right of the cutting line. 

In the high energy limit s — > oo, A — > and T,^\b) = T A (b). It is easy to see 
that for j > 1 the last two terms in the square bracket of cr^') cancel with the first 
two terms of We are thus left with the first two terms of the last two 

terms of and a£\ This gives 

o?(b) = of (a + a)- at -5) + a(l - e)) . (9) 

Since 5 = 1 the last term does not contain any collision of the projectile. Only a part 
of the first term corresponds to the J/ip production mechanism we have assumed. 
More precisely, we have to consider only the inelastic part of o^ia + a). Moreover, 
as in eqs. (2) and (3), we have to remove the terms in which only a and/or a N °^ 
are present. We thus have 



k=l 



°» = e(;) (*+*) k -{<f+t) 



n{b) 



[l-(a + a)T A (b)] A - k = ~ Aat N T A (b) (10) 

which reproduces the result in eq. (3). 

Let us now turn to the low energy limit in which the coherence length tends to 
zero (A — > oo). In this case the only surviving contribution is cr^ =1 ) - the only one 
where the cosinus damping factor is absent. Note that T^ =1 \b) = T A {b). Removing 
the second term (as) in the first factor of (8) (which corresponds to no interaction 
of the projectile since 8 — 1) we get in this limit 

a pA (0) = <?pN — — ( 1 1 ) 

(Tabs 

with <j a i> s — a(l — e) + a(l — 5). Since 8 — 1 we obtain the conventional Glauber 
formula with a rescattering controlled by a ahs = a(l — e) - which is present for J/ip 
production (e ~ 0) and absent for open charm (e = 1). 

So far our formulae are valid for mid-rapidities. If we want to consider large 
positive values of xp, we have to take into account that the geometry introduces a 
correlation between xf and impact parameter b. Indeed, for very peripheral collision 
(where the projectile goes through the edge of the nucleons) the number of wounded 



nucleons (i.e. cut interactions of the projectile with target nucleus) is small and 
production at large xp is favored. In ref. [7] this correlation was taken into account 
phenomenologically by introducing a softening factor = (1 — xj) k ^ 1 . Here x\ (x 2 ) 
is the x variable relative to the projectile (target) with x\x 2 s = M 2 Z and xf = x\—x 2 - 
k is the number of cut interactions (in Ref. [7], 7 = 2 for J ftp production). With 
the factor F k the spectrum of J/ip gets softer with increasing k. Since in eqs. 
(7), (9) and (11) we have kept track of the cut interactions it is immediate to 
introduce the softening factor F k . In the low energy limit we have to replace a (5—1) 
by a 5(1 — xj) — o = —ox\ (5 = 1). In this way, the absorptive cross-section 
(fabs = (1 — e)a is replaced by an effective cross-section a e ff = a abs + ax\. Thus, 
when x\ — > 1 we get rescatterings controlled by a typical hadronic cross-section. It 
was shown in [7] that such a large cross-section is needed in order to reproduce the 
J/ip suppression at large x F . Note that the same increase of the effective absorptive 
cross-section as x\ — > 1 occurs for open charm production. Likewise, the softening 
factor has to be introduced in the high energy limit. However, its effect at the highest 
available rapidity at RHIC energies (y = 2.2) is very small and can be neglected. 

The extension of the above formalism to AB collisions is straightforward (see for 
instance Ref. [8]). The formulae are obtained by a convolution of the ones for pA 
and pB. 

3 Nuclear shadowing 

The above considerations -that is to say, the A 1 dependence in eqs. (3) and (10)- 
suggest a phenomenological analysis of J/ip suppression in dAu collisions at RHIC 
energies using a a bs = and taking only into account nuclear shadowing. The latter is 
particularly important because the data are presented as a ratio <J J Jau/N C oII ~ rather 
than a ratio ^' J JauI a dAu m "which shadowing effects would cancel to a large extent. 

We call shadowing the mechanisms that makes that the nuclear structure func- 
tions in nuclei are different from the superposition of those of their constituents 
nucleons. Typically, the shadowing behaves as follows: it increases with decreasing 
x and decreases with increasing Q 2 . Several explanations to the shadowing have 
been proposed. We can distinguish between two main approaches: models based on 
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multiple scattering [8]- [10] in the rest frame of the nucleus, and models based on a 
frame in which the nucleus is moving fast, and where gluon recombination due to 
the overlap of the gluon clouds from different nucleons reduces the gluon density in 
the nucleus. 

In the rest frame on the nucleus, nuclear shadowing can be seen as a consequence 
of multiple scattering [8]- [10] the incoming virtual photon splits into a colorless qq 
pair long before reaching the nucleus, and this dipole interacts with typical hadronic 
cross sections which results in absorption. Multiple scattering can be related to 
diffraction by means of the Abramovsky-Gribov-Kancheli (AGK) rules. 

On the other hand, another possible approach consists on the models based on 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution [11] like EKS98 [12], 
nDS [13] or HKN [14]. This type of models do not try to address the origin of nuclear 
shadowing but study the (Revolution of nuclear ratios of parton densities, 

R ^Q 2 ) = A0^gTy fi = 9>9>9, (12) 

through the DGLAP evolution equations. They try to perform for the nuclear case 
the same program developed for the nucleon: nuclear ratios are parameterized at 
some value Ql ~ 1^-2 GeV 2 which is assumed large enough for perturbative DGLAP 
evolution to be applied reliably. These initial parameterizations for every parton 
density have to cover the full x range < x < 1. Then these initial conditions are 
evolved through the DGLAP equations towards larger values of Q 2 and compared 
with experimental data. From this comparison the initial parameterizations are 
adjusted. While DGLAP approaches do not address the fundamental problem of 
the origin of shadowing, they are of great practical interest. They provide the parton 
densities required to compute cross sections for observables characterized by a hard 
scale for which collinear factorization can be applied, 

In the present paper, we are going to compare the data with the results obtained 
from [8] -model based on multiple scattering- and from EKS98 [12] -model based on 
the GRV LO parton densities and where the shadowing ratios for each parton type 
are evolved to LO for 1.5 < Q < 100 GeV-. 

Following the first approach, in order to compute the shadowing we need the to- 
tal contribution which arises from cutting the two-exchange amplitude in all possible 
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ways (between the amplitudes and the amplitudes themselves in all possible man- 
ners). It can be shown that, for purely imaginary amplitudes, this total contribution 
is identical to minus the contribution from the diffractive cut. Thus diffraction be- 
comes linked to the first contribution of the nuclear shadowing. 

In consequence, we apply here a dynamical, non linear shadowing [8] which is a 
consequence of multiple scattering, and it is determined in terms of the diffractive 
cross-sections. It would lead to saturation at s — > oo. It is controlled by triple 
Pomeron diagrams. This kind of shadowing gives a positive contribution to the 
diffractive cross-section, and, on the contrary, its contribution to the total cross- 
section is negative. It gives good results for the kinematic regions where nuclear 
DIS data exist [15]. The effect of the shadowing corrections in our approach [8] -[9] 
can be expressed by the factor: 

n AB(b) = l£±MLMz£> (13) 

J-abkP) 

where 

m») - r^mzm ■ (14) 

Here the function F(s) is given by the integral of the ratio of the triple Pomeron 
cross-section d 2 a ppp / dydt at t — to the single Pomeron exchange cross-section 
cr p (s) : 

* 



rDmax 

F(s,y)=4ir dy — 

J y rain ® p 



1 d 2 ~ PPP 



<T p (s) dy dt 



= C [exp (Ay max ) - exp (Ay min )} (15) 

t=o 

with y = £n(s/M 2 ), where M 2 is the squared mass of the diffractive system. For 
a particle produced at y cm = 0, y max = \tn(s/m T ) 2 and y min = £n(R A m N /V3). 
A = op(0) — 1 = 0.13 and C is a constant proportional to the triple Pomeron 
coupling. Ra is the nuclear radius, T^(6) the nuclear profile function and Tas(6) = 
/ d 2 sTA{s)TB{b—s). The denominator in Eq. (13) correspond to the sum of all "fan" 
diagrams with Pomeron branchings (generalized Schwimmer model [16]). Note that 
shadowing corrections to inclusive spectra are not specific to soft processes. The 
triple Pomeron terms described above are also responsible for shadowing in hard 
processes. 

We have then for the suppression of hadron h production in AB collisions 

Sk{b > S ' V) = 1 + AF h (y A )T A (s) l + BF h (y B )T B (b-s) (16) 
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where F h (y A ) = C[exp(Ay max ) -exp(Ay mi „)], y min = £n(R A m N / A = 0.13 and 
C = 0.31 fm 2 . y max depends on the rapidity of h : y^ ax = \£n{s/m^) ± y with the 
+ (— ) sign if h is produced in the hemisphere of nucleus B(A). itlt is the transverse 
mass of h. For charged particles we use m T = 0.4 GeV and for J ftp ttit = 3.1 GeV. 
For dAu collisions we just take the first factor of (16), which corresponds to pA 
collisions. With a abs = this suppression equals the total J/ip suppression, since no 
hot -final state interactions- effects are expected in this case. 

We are going to compare the results of our model with the ones obtained from the 
EKS98 parameterization. In order to do that, we have applied the same procedure 
as the one introduced by R. Vogt in Ref. [17]. The spatial dependence can be 
parameterized assuming that shadowing is proportional to the local density, pa(s), 

R?(x,Q 2 ,r,z) = l + N ws [R?(x,Q 2 ) - , (17) 

Po 

where r and z are the transverse and longitudinal location in position space, with 
s = \J\r\ 2 + z 2 , Pa( s ) corresponds to the Woods-Saxon distribution for the nu- 
cleon density in the nucleus, po is the central density, given by the normalization 
/ d 2 rdzp A (s) = A and Rf(x,Q 2 ) is the shadowing function from EKS98 as defined 
by eq. (12). N ws is chosen so that (I /A) J d 2 r J dzp A (s)Rf(x,Q 2 ,f,z) = Rf(x,Q 2 ) 
In the small x limit, the above formulae transforms into 

R?(x, Q 2 , f, z) = l + N p [Rf(x, Q 2 ) - l] { d f PA ^ Z \ ■ (18) 

J dzp A (U,z) 

The integral over z includes the material traversed by the incident nucleon, so we are 
considering that the incident parton interacts coherently with all the target partons 
along its pathlength. 

The normalization requires again (l/A) J d 2 r J dzp A (s)Rf(x,Q 2 ,r,z) = Rf(x,Q 2 ). 
We do not have included the nuclear absorption, since we consider a ahs = 0. 



4 Results 

Our results of the nuclear modification factor for the J/ip -that is, the AB/pp ratio 
for the JM, dN jfJ dy , at a given centrality, normalized by the number of collisions 

dN p y /dy 

N^ou at this centrality, R A b = N AB^J/^ d — ^ or dAu collisions, compared with 
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PHENIX data [2], are presented in Fig. 1. They are also compared with the results 
obtained in the framework of the EKS98 model [12, 17]. 

The agreement between theory and experiment is quite good in both cases. Note 
that the increase of J / ip suppression with rapidity is entirely due to the correspond- 
ing increase of shadowing. 

The situation is different in Au Au collisions. Indeed, it can be checked with our 
formulae that shadowing effects are, in this case, practically independent of rapidity. 
This is due to the opposite signs of y in y^ ax and y^ iax - As a consequence, the J/ip 
suppression in Au Au in a comovers approach - a charmonium state produced in a 
primary nucleon-nucleon collision is dissociated through interactions with the dense 
medium subsequently formed in the collision - is maximal at y — 0, where the 
comovers density is maximal, and decreases when we move to positive or negative 
rapidities. 

Moreover, this is the case in any model that considers the suppression of the 
J/ijj dependent on the medium density, that is, such dissociation could occur in a 
confined [18] as well as in a deconfined medium [19]. 

One possibility to reverse this tendency should be the recombination of c-c pairs 
into J/ip [20]. In this case, in the hadronization of the quark-gluon plasma, charmo- 
nium formation can occur by binding of a c with a c from different nucleon-nucleon 
collisions, as well as from the same. If the total number of available c-c pairs consid- 
erably exceeds their thermal abundance, statistical recombination enhances hidden 
relative to open charm production, as compared to hadron-hadron collisions. This 
possibility has been neglected in our model. It would correspond to the introduction 
of a gain term in differential equations which govern the final state interactions. We 
will develop this possibility in our future work. 

The results for Au Au collisions at RHIC at y = with cr a b s = and o co = 
0.65 mb (the same comovers interaction cross-section obtained from the CERN data) 
were given in [8]. We present them again in Fig. 2 and compare them with the 
PHENIX data [21]. The corresponding results at y — 1.7 are also presented. The 
agreement with the data is considerably improved by taking a a b s = 0. Nevertheless, 
our results are still somewhat lower than the data at y — 0, while there are higher at 
y = 1.7. This shows that in the present version of our model the rapidity dependence 
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of the suppression is not reproduced. Moreover, it is in contradiction with the present 
data, as it has been indicated above: the J/ip suppression in Au Au is maximal at 
y — 0, where the medium density is maximal, and decreases when we move to 
positive or negative rapidities. 

Note that a co is the comover cross-section properly averaged over the momenta 
of the colliding particles and over the collision time. A large contribution to this 
interaction comes from the few first fm/c after the collision - where the system is 
in a pre-hadronic or partonic stage. Actually, Brodsky and Mueller [22] introduced 
the comover interaction as a coalescence phenomenon at the partonic level. In view 
of that, there is no precise connection between a co and the physical J/ip — n or 
J/ip — N cross-section, and a co has to be considered as a free parameter. Since it 
is a cross-section near threshold, we do not expect a substantial energy variation. 
Because of this, we have used the same comovers interaction cross-section obtained 
from the CERN data. However, this cross-section does not have to be identical at 
CERN and RHIC - since the momentum distribution of the comovers can change 
with energy. Unfortunately we are unable to evaluate such an effect. Another source 
of uncertainty in our calculation resides in the comovers density. As an illustration, 
we show in Fig. 2 the effect on the J/ip suppression at y = resulting from a 20% 
decrease in this density. 

More important, theoretical calculations [23] show that the energy dependence of 
the J / ip-h&droia cross-section increases very fast with increasing energy near thresh- 
old. Because of that, the decrease of the comovers density at forward rapidities 
might be (over) compensated by an increase of the comovers cross-section a co . In 
this case, the maximum of the J/ip suppression may take place at a value of the 
rapidity different from y = where the comovers density is maximal. This point, as 
well as the possibility of recombination, are at present under investigation 3 . 

5 Conclusions 

In conclusion, we have presented theoretical arguments supporting the idea that, 
at high energy, nuclear absorption of the J/tp vanishes (dabs = 0). Only standard 

3 We thank A. B. Kaidalov for interesting suggestions on this point. 
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nuclear effects (shadowing) are present. This produces an increase of the J ftp sup- 
pression with increasing rapidity in dAu collisions. The results at \fs = 200 GeV 
energies can be described in this way. In Au Au collisions the shadowing effect is 
practically independent of rapidity and the J/ip suppression in the comovers inter- 
action model presented here decreases when moving away from mid-rapidity, due to 
the corresponding decrease of the comovers density. 
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Figure 1: (Color online) Results on suppression for dAu collisions vs. the 
number of collisions in 3 different rapidity regions. We take into account nuclear 
shadowing alone (<r a bs — 0). Data are from PHENIX. Continuous line: our results 
from Pomeron shadowing, discontinuous line: results from EKS98 shadowing. 
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Figure 2: Results on J/tp suppression for AuAu collisions vs. the number of par- 
ticipants at y = (solid line) and y — 1.7 (dotted line). The dashed-dotted line is 
the result at y — obtained by decreasing the comovers density by 20% (see main 
text). The dashed lines correspond to the suppression from shadowing alone. Data 
are from PHENIX (black points: data at y — 1.7, open points: data at y — 0). 
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